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Abstract 

We apply the one-loop results of the SU(3) L x SU(3) R ChPT su- 
plemented with the inverse amplitude method to fit the available ex- 
perimental data on tttt and irK scattering. With esentially only three 
parameters we describe accurately data corresponding to six differ- 
ent channels, namely (I, J) = (0, 0), (2, 0), (1, 1), (1/2, 0), (3/2, 0) and 
(1/2,1). In addition we reproduce the first resonances of the (1,1) 
and (1/2, 1) channel with the right mass corresponding to the p and 
the K*(892) particles. 
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Introduction: In 1979 Weinberg [1] suggested that it is possible to 
summarize many previous current algebra results in a phenomenological la- 
grangian that incorporates all the constraints coming from the chiral sym- 
metry of the strong interactions and QCD. This technique, also called Chiral 
Perturbation Theory, was developped some time later to one-loop level and 
in great detail in a celebrated set of papers by Gasser and Leutwyler [2]. In 
these works, the authors showed how it is possible to compute many differ- 
ent Green functions involving low-energy pions and kaons as functions of the 
lowest powers of their momenta, their masses and a few phenomenological 
parameters. In these and in further works it was also shown that this method 
provides a good parametrization of many low-energy experimental data. 

More recently, the strongly interacting symmetry breaking sector of the 
Standard Model [3] has also been described phenomenologically using ChPT 
[4]- 

In both contexts , one of the main obstacles found when one tries to ap- 
ply ChPT to higher energies relies in the issue of unitarity. ChPT, being 
a consistent theory, is unitary in the perturbative sense. However, as the 
expansion parameters are the momenta and the Nambu-Goldstone boson 
masses, perturbative unitarity breaks down, sometimes even at moderate en- 
ergies [5]. Different attempts to improve this behaviour of ChPT and extend 
the applicability to higher energies, have been proposed in the literature. 
These methods include the use of unitarization procedures such as the Pade 
expansion [5,6], the inverse amplitude method [5], and the explicit introduc- 
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tion of more fields describing resonances [7]. All them improve the unitarity 
behavior of the ChPT expansion and provide a more accurate description of 
the data although there is some controversy about which of them is more 
appropriate. 

In this note we show the results of the application of the inverse amplitude 
method to the ChPT one-loop computation of elastic irir and ttK scattering 
[2,8]. The obtained amplitude converges at low energies with that of stan- 
dard ChPT but it fulfills strictly elastic unitarity, since the contribution of 
the right cut in the corresponding dispersion relation is taken into account 
exactly in the inverse amplitude method. Using this approach (which inci- 
dentaly is connected with the formal [1,1] Pade approximant of the one-loop 
amplitudes) we will make a global fit of the data for elastic ttti and irK scat- 
tering just varying the L ± , L 2 and L 3 ChPT constants. The resulting values 
of these parameters will not be quite different from those obtained previ- 
ously, but the range of energies and the quality of the fit will be enlarged 
quite amazingly. 

The partial waves in ChPT: The elastic scattering partial waves are 
defined from the corresponding isospin amplitude Tj(s,t) as: 



where K — 2 or 1 depending on whether the two particles in the reaction 
are identical or not. For elastic tttv scattering the possible isospin channels 
are / = 0, 1, 2 whilst for irK we can have I = 1/2 or J = 3/2. In the first 
case the isospin amplitudes T/ can be written in terms of a simple function 
A(s, t, u) as follows; T (s, t, u) = 3A(s, t, u)+A(t, s, u)+A(u, t, s), T^s, t, u) = 




(1) 
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A(s, t, u) — A(t, s, u) and T 2 (s, t, u) = A(s, t, u) + A(t, s, u). In the second, 
the 1=1/2 amplitude can be written as T 1 / 2 (s,t,u) = (3/2)T 3 / 2 (m, t, s) — 
(l/2)T 3/2 (s, t, u). The Mandelstam variables s, t, u fulfill s + 1 + u = 2{M 2 a + 
Mj) where we use the notation ct and (3 = n or K, so that we can describe 
both processes with the same general formulae. For s > s t h = (M a + Mp) 2 
the partial waves tu can be parametrized as: 

t IJ (s) = e«"WsHi6 I j(s)/<T a f ) (s) (2) 

where 

<W0 = y(l " — " ) (3) 

The tu(s) amplitude fulfills the elastic unitarity condition 

Imtu = a a(3 \t u \ 2 (4) 

on the right cut. 

Using standard one-loop ChPT it is possible to compute the above scat- 
tering amplitudes to order p A (as it is customary, any Mandelstam variable, 
Ml and M\ will be considered of the order of p 2 ). The relevant functions 
A(s,t,u) and T 3 / 2 were computed in [2,8]. 

A(s,t,u) = (s-M 2 )/F 2 + B(s,t,u) + C(s,t,u) + 0(E 6 ) (5) 

1 M 4 1 1 

B(s,t,u) = T -S^J^s) + -{s 2 -Mt)JlAs) + -/J r KK{s) 

+ ±(f - 2Ml) 2 J^{t) + t(s - u)[M^{t) + \w KK {t)\ + (t~ «)) 

C(s,t,u) = ^mi + L 3 )(s-2M 2 ) 2 + U 2 [(t~2M 2 ) 2 + (u-2M 2 ) 2 ] + 

TV 

+ (ALl + 2Ll)M%s - 2M 2 ) + (8L r 6 + 4L r 8 )M*) 
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and 

M 2 4- /If 2 — <? 

T^(s,t,u) = - + 2F f + Tl(s,t,u) + T[(s,t,u) + T^(s,t,u) (6) 



+ ^A-[2A^ + 6M^-4s]+^[5Ag + 7M^-6s]) 
2 

"2 



+ 2U 2 [{s-Ml-Mlf + {u-Ml-Mlf} 

+ ^[(w-M2-M|) 2 + (t-2M2)(t-2M|)] 

+ 4L^(M* + M 2 K ) - 4M*M£] 

+ 2LlM 2 w (M 2 w -M 2 K -s) + 8(2L r 6 + Ll)M 2 w M 2 K ) 

Tf(s,t,u) = ^^(t(u-s)[2M r wn (t) + M r KK (t)} 

+ 1[( 8 - t)(L wK (u) + L Kv (u) - u(M r nK (u) + M r Krj (u))) 

+ {M 2 K - M 2 )(M: k (u) + M r Kv (u))} + J: k (s)(s -Ml- Ml) 2 

+ 1 -{M 2 K - M 2 n )[K nK (u)(5u - 2M 2 K - 2M 2 ) + K Kv (u)(3u - 2M 2 K - 2M 2 )} 

+ l -Jl K {u)[llu 2 - Y2u{M 2 K + Ml) + 4(M| + Ml) 2 ] 

+ \Jk,{u){u - \{M 2 K + Ml)) 2 + Ij^tm ~ Ml) 

+ \jWit)t 2 + \r m {t)Ml{t- % -M 2 K )) 



The masses M a and the decay constants F a appearing in these equations are 
the physical values. The relation with the corresponding constants appear- 
ing in the chiral lagrangian, and the functions fi a can be found in [8]. The 



transcendental functions M*p, L r a/3 and are defined in [2]. The first terms 
in the above amplitudes reproduce the well known Weinberg low-energy the- 
orems. The L\ constants can be considered as phenomenological parameters 
that up to constant factors are the renormalized coupling constants of the 
chiral lagrangian renormalized conventionaly at the m v scale. Their rela- 
tion with the corresponding bare constants Lj and their evolution with the 
renormalization scale can be found in [2]. Using eq.l and eq.5 it is possible 
to obtain the corresponding partial wave amplitudes. In the general frame- 
work of ChPT they can be obtained as a series with increasing number of p 2 
powers, i.e. 

tu = t<?} +$ + ■■■ (7) 
where tfj is of order p 2 and corresponds to the low-energy theorem and tf] 
is of order p 4 . In general, the real part of tfj cannot be expressed in terms of 
elementary functions but it can be computed numerically. The amplitudes 
in eq.5 have been used in the literature, without further elaboration, to fit 
the low-energy tttt and ttK scattering data [2,8,9]. 

Dispersion relations and the inverse amplitude method: One 
very important point concerning the partial wave amplitudes computed from 
ChPT to one-loop is the fact that they have the appropriate cut structure, 
namely the left cut and the right or unitarity cut. However, they only fulfill 
the unitarity condition on the right cut in a perturbative sense i.e. 

Imtfj = a aP | tfj | 2 (8) 

Let us show now how, with the use of dispersion theory, it is possible to 
build up a completely unitarized amplitude for the ttti and irK amplitudes 



starting from the one-loop ChPT result above. Let us start writing a three 
subtracted dispersion relation for the partial wave tij\ 

^ ^ ^ 9 s 3 r 00 Imtrj(s')ds' 
tu(s) = C + C lS + C 2 s 2 + - / - IJK ' , + LC{t LJ ) 9 

7T J{M tx +M g ) 2 S' A {S' — S — 16) 

where LC(tjj) represent the left cut contribution. In particular, this equation 
is fulfilled by the standard one-loop ChPT result. Note that three subtrac- 
tions are needed to ensure the convergence of the integrals in this approx- 
imation since the one-loop ChPT amplitudes are second order polynomials 
modulo log factors (higher order ChPT amplitudes require more subtractions 
but we do not know how many are needed for the exact result). To the one- 
loop level, the imaginary part of the amplitude can be written on the right 
cut integral as Imtu ~ Imtf] = utfj 2 (from the left cut there is a generic 
tf) contribution). 

The substraction terms appearing in the above dispersion relation can be 
expanded in powers of M^, and they have contributions from tfj and from 
tf] that can be written as a + ais and 6 + bi s + ^s 2 and also from higher 
order terms of the ChPT series i.e Cq = ao + b + ... and C\ = ai + b± + .... The 
a and b constants depend also on the L\ parameters. Then, in the one-loop 
ChPT aproximation, the exact amplitude t(s) can be written as: 

t u ( S )^t?]( S ) + t?)(s) (10) 

where: 

4j — a o + a i s (11) 
4j = b + b lS + b 2 s 2 + 



7T J(M a +M l3 ) 2 S ,a {s' — S — 16) 



where in the left cut contribution we have to integrate here Imt I j{s) as 
an aproximation to tu(s). In some sense, we can understand eq.9,10 and 
eq.ll in a rather different way that will be useful later: we can assume that 
the dispersion relation in eq.9 is fulfilled by the exact amplitude t(s) and we 
solve this equation approximately by introducing inside the left and right cut 
integrals the one-loop prediction for Imt(s) to find again eq.10. However, we 
would like to stress again that this result does not fulfill the elastic unitarity 
condition in eq.4 but only the perturbative version in eq.8. In fact this 
happens to any order in the ChPT expansion since a polynomial can never 
fulfill eq.4. However, there are other ways to use the information contained in 
the ChPT series apart from the direct comparison with the experiment of the 
truncated series. Instead of using the dispersion relation for the loop ChPT 
amplitude we can try a three subtracted dispersion relation for the inverse 
of tjj or more exactly for the auxiliar function G(s) = tfj 2 /tu namely: 

G(s) = G + G 1 3 + G 2 s 2 + - r l™ G{s ' )ds [ +LC{G)+PC (12) 

7T J(M ct +M /3 ) 2 S' S (S' - S - 16) 

where LC (G) is the left cut contribution and PC is the pole contribution 
that eventually could appear due to possible zeros of tjj(s). Now, on the 
right cut we have ImG = tfj 2 Im(l/tjj) = —tfj 2 Imtjj/ \ tjj | 2 = —tf) 2 a. 
This means that the right cut integral appearing in the dispersion relation 
for G is the same than that appearing in eq.10 and eq.ll for the one-loop 
dispersion relation for t and hence it can be obtained from the one-loop 
result. The left cut integral of the dispersion relation in eq.12 cannot be 
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computed exactly but we can use the one-loop ChPT result to write ImG = 
—tfj 2 Imtu/ | tu | 2 ~ —Imtfj. Besides, the substraction constants can be 
expanded in terms of M^/Fj powers so that G = a - b + 0((M 2 /F|) 3 ), 
F^G 1 = a 1 -6 1 + 0((M2/F|) 2 ), and F*F*G 2 = -b 2 + 0{M 2 J Fj). Therefore, 
neglecting the the pole contribution, the dispersion relation in eq.12 can be 
written approximately as: 



t 



(0)2 

~ a + ai-s — b — bi-s — b 2 s 2 (13) 
f°° atf} 2 (s')ds' 

{Mc+Mp) 2 S ,3 (s' — S 



vr J(M a +Mf,) 2 s' A [s' — s — ie) 



with LC(G) computed with the ImG approximated by —Imtfj or in other 
words: 

(0)2 

v 



or what it is the same: 



(14) 

^(0)2 

* WZTw (15) 

Remarkably, to derive this result the one-loop ChPT approximation has 
been used only inside the left cut integral but not inside the right cut integral 
which was computed exactly. This is in contrast with the one-loop ChPT re- 
sult in eq.10 when considered from the point of view of the dispersion relation 
that tjj has to fulfill. In this case, in adddition to the above approximations 
used to derive eq.15, we have done a much stronger one concerning the right 
cut contribution since the Imtu was also approximated by the one-loop re- 
sult. This fact is crucial because the right cut contribution is responsible 
for the very strong rescattering effects present in these reactions. As a con- 
sequence of that, eq.15 fulfills the unitarity condition Imtjj = a \ tjj 
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exactly and not only perturbatively as it was the case of eq.10. Of course 
both results (eq.10 and eq.15) provide the same answers at low energies but 
it is expected that eq.15 will provide more realistic results at higher energies 
and this is in fact what we find when comparing with the experimental data. 

When zeros are present in the partial waves, they need to be included 
and the final result is not so simple as in eq.15. This can be done taking into 
account their contribution to eq.12 or just making a new sustraction in these 
points. However, in practical cases, it can happen that the corresponding 
residuous are very small and eq.15 is still a good approximation outside the 
region around the position of the these zeros. This is for instance the case 
of the I = J = and I — 2, J — channels [6] which have zeros close to 
and 2Mfj- respectively. However, for a given channel I J there is no way 
to know a priori if zeros will be present or not since ChPT only provides a 
low-energy expansion of t(s) and not the whole amplitude. Nevertheless, in 
practice, a simple inspection of the low-energy behaviour of the amplitudes 
can make sensible the hipothesis that zeros are not present in the low-energy 
region but of course this must not be necessarily true for all the channels. 

We could also ask about the possibility of extending the above method to 
higher orders of the ChPT series for instace to the two-loop computation. In 
principle this can be done in a straightforward way. We start writing a four 
subtracted dispersion relation for the two-loop amplitude since it is a third 
order polynomial modulo logarithms. As in the one-loop case, this dispersion 
relation can be interpreted as an approximation to the exact amplitude tjj. 
Then we write a four subtracted dispersion relation for the auxiliar function 
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G = tjj jtij. The integrals of this dispersion relation can be related modulo 
higher order terms to the ChPT by approaching ImG by its low energy 
expansion. The subtraction constants can be treated similarly as we did in 
the one-loop case an finaly we arrive to: 

t (0)2 

tlJ ~ W) M) 7 I (1)2 /f (o) 7m ( 16 ) 
l IJ L IJ L IJ I l IJ L IJ 

It is very easy to show that this amplitude, which incidentally corresponds 
to the formal [1,2] Pade approximant, fulfills the elastic unitarity condition 
in eq.4 but in this case, the right cut integral appearing in the dispersion 
relation for G cannot be computed exactly as it was in the one-loop case. 
Unfortunately, as there is not any two-loop computation available, we cannot 
confront the above equation with the experimental data. 

Determination of the chiral parameters and discussion: First of 
all, we have used the set of low energy chiral parameters Lj proposed in [11] 
for i=l,2,3 and those given in [2] for i=4..8 which will not be changed in 
this work. Appliying those values to the proposed unitarized partial waves 
of eq.15, we have found the appropriate qualitative behaviour for both re- 
actions. Moreover, we find two resonances in the partial wave amplitudes 
corresponding to the (1,1) and (1/2,1) channels. These resonances can be 
identified naturally as the p and the K* particles. The phase-shifts cross the 
7r/2 value in their corresponding channels. In addition, in the chiral limit 
where an analytical computation of the partial waves is possible, two poles 
are found in he second Riemann sheet accordingly to these resonances. In the 
other channels, where no physical resonances exist below lGeV, they do not 
appear in our numerical results. We consider this fact as a strong support 
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for the use of the unitarized partial waves in eq.15. 

To fit completely the data with these formulae the next step has been 
to tune the (1,1) irn channel to give the correct mass for the p resonance. 
Since this partial wave is almost only sensible to the relation L 3 + 2L\ — L 2 , 
then, setting the p mass to 774 Mev means fixing a value of this special 
combination of parameters. 

Nevertheless, there are still two degrees of freedom. Slightly varying 
the initial values we can fit the nn (2,0) and (0,0) channels, thus obtaining 
L\ and L 2 . Unfortunately the experimental data coming from these chan- 
nels allow some small uncertainty in the parameters, rending a K* mass 
between 850 and 950 Mev. So we use the mass of this resonance for a fur- 
ther parameter tuning, and finally we obtain, at the renormalization scale 
p — M v — 548.8Mew: 

L\ = 0.610" 3 , U 2 = 1.610 -3 , L 3 = -3.810 -3 (17) 

These values are well inside the errors quoted in [11] 

L\ = (0.88±0.47)10' 3 ,L^ = (1.61±0.38)10~ 3 , L 3 = (-3.62±1.31)10 -3 (18) 

The K* mass thus obtained is 880 Mev In figures 1 and 2 can be seen 
the results of this global fit for tttt scattering and those from non unitarized 
ChPT with the U parameters proposed in [10]. Figures 3 to 5 represent our 
global fit curves (continuous lines) for nK scattering in three different chan- 
nels. The dashed lines are the non unitarized ChPT predictions from the Lj 
given in [2] and [11]. Fig. 3 is the (3/2,0) channel and figures 4 and 5 are the 
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(1/2,1) (1/2,0) respectively. 



Conclusions: The inverse amplitude method applied to the one-loop 
result coming from ChPT produces a simple way to unitarize the Goldstone 
boson elastic scattering amplitudes which takes into account, exactly, the 
strong rescattering effects. Incidentally, this method is formally equivalent 
to the [1, 1] Pade approximant applied to the one- loop ChPT result, pro- 
vided that the exact partial wave amplitude has no zeros in the first sheet. 
Note that since the one-loop result is not, strictly speaking, a polynomial, the 
equivalence is only formal (This is not the case if one considers the amplitudes 
as polynomials in i ) 

The unitarized amplitudes (with the previously fitted parameters for the 
standard one-loop ChPT result) give rise to the appearence of two resonances 
in the (1,1) and (1/2,1) channels that have to be understood as the p and 
K*, but not any more are found in other channels where no physical reso- 
nances exist below lGeV. Therefore, the existence of these resonances is a 
highly non trivial prediction of the approach followed here and do not need 
to be introduced by hand in the data fit . Just tunnig slightly only three 
parameters, namely L 1? L 2 and L 3 we are able to obtain the right value for 
the masses of the p and K* resonances. In addition we provide a fit for six 
channels in a remarkable agreement with the experiment. 

The results obtained in this and previous works strongly suggest that 
the range of validity of the one-loop ChPT can be enlarged by a consistent 
traitment of the analiticity and unitarity constraints summarized in the dis- 
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persion relations but a great deal of work still remains to be done in this 
direction. As a final comment, we think that the recently proposed large 
N approximation to ChPT (with N being the number of Nambu-Godstone 
bosons) [25] (see also [26]) could provide a new extension of the one- loop 
ChPT results (see also [26] for applications to the 77 — > nn reaction [27]). 
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Figure Captions 

Figure 1 .- (1,1) Phase shift for irir scattering. The continuous line corre- 
sponds to our fit using eq.15. The dashed line is the result coming from non 
unitarized ChPT with the l{ parameters proposed in [10]. The experimental 
data comes from: O ref. [16], A ref.[20]. 

Figure 2 .- Phase shift for tt-k scattering. The results coming from the 
fit proposed in this paper (eq.15) are the continuous line which represents 
the (0,0) phase shift, and the dashed line which corresponds to that of (2,0). 
The dotted and dashed-dotted lines are the (0,0) and (2,0) phase shifts re- 
spectively,they were obtained with non unitarized ChPT and the parameters 
given in [10]. The experimental data corresponds to: A ref.[12], O ref.[13],D 
ref.[14],0 ref.[15], y ref.[16],* ref.[17],x ref.[18],« ref.[19]. 

Figure 3.- Phase shift of the (3/2,0) channel for irK scattering. The 
continuous line is the result of the inverse amplitud method (eq.15) with the 
parameters proposed in this paper, whereas the dashed line is non unitarized 
ChPT with the parameters proposed in [2] and [11]. Data corresponds to 
ref.[21] 

Figure 4.- The same as figure 3 but for the (1/2,1) partial wave. Data 
comes from: A ref. [21] and * ref.[24]. 

Figure 5.- As figure 3 but for the (1/2,0) channel. The experimental data 
corresponds to: A ref.[21],<) ref. [22], and * to ref. [23]. 
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